Steady state entanglement of two atoms created by classical driving field 
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The stabilization of steady state entanglement caused by action of a classical driving field in 
the system of two-level atoms with the dipole interaction accompanied by spontaneous emission is 
discussed. An exact solution shows that the maximum amount of concurrence that can be achieved 
in Lamb-Dicke limit is 0.43, which corresponds to the entanglement £max = 0.285 ebit. Dependence 
of entanglement on interatomic distance and classical driving field is examined numerically. 
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I. INTRODUCTION 

It is well known that quantum entanglement plays fun- 
damental role in the quantum information processing and 
quantum computing [1|, ^ . The practical application of 
entanglement requires the robust entangled states. This 
notion includes long enough lifetime of the states and 
high amount of entanglement (as close to maximum en- 
tanglement as possible). 

The simplest entangled states of two or more qubits 
can be modelled physically by the states of two-level 
atoms, which have been successfully used for decades as 
the main tool for testing fundamentals of quantum me- 
chanics 3]. In particular, a number of important exper- 
iments on generation and manipulation with entangled 
states in atomic systems have been performed recently 
(e.g., see 4] and references therein). 

The lifetime of entangled states in atomic systems is 
usually governed by the spontaneous emission rate and 
therefore is quite short Q. The problem of generation 
of robust entangled states in atomic systems is of high 
importance. In fact, this problem belongs to a more gen- 
eral class of problems related to the so-called quantum 
decoherence. This branch of modern physics examines 
an irreversible evolution of quantum systems under influ- 
ence of a dissipative environment (see la| and references 
therein). Thus, using the ideas of quantum decoherence 
theory, we have to choose a proper environment for an 
atomic system, providing an irreversible evolution of this 
system towards a stable (or at least metastable) entan- 
gled state. 

For example, it has been proposed in Refs. |3,la| to use 
the three-level atoms with A-type configuration of levels 
instead of the two-level atoms. In this case, an unstable 
maximum entangled state is created by absorption of a 
photon, resonant with the transition between the ground 
and highest excited level in the system of the two atoms. 
This state rapidly decays either to the ground state or to 
the intermediate excited state. The latter process is ac- 
companied by generation of the so-called Stokes photon 
and leads to creation of a new maximum entangled state 
of two atoms. If the Stokes photon, generated by the 
transition between the highest and intermediate excited 
states, is then discarded through the use of an appropri- 
ate dissipative environment, the system remains in the 



maximum entangled state for a long time determined by 
the non-radiative decay of the intermediate atomic state. 

In the above example, the maximum entangled atomic 
state is induced by the photon exchange between the 
atoms. Another physical proposal is based on the use 
of direct dipole interaction between the two-level atoms 
0. In this case, the energy loss in the system due to the 
spontaneous emission is supposed to be compensate by 
interaction with an environment. In Refs. Q it was pro- 
posed to use the interaction of atoms with the squeezed 
vacuum state in addition to the dipole-dipole interaction. 
It was shown that a steady state with reasonable amount 
of entanglement can be achieved in this case. 

This is not an unexpected result. The point is that it 
has been recognized recently that entanglement is pecu- 
liarly connected with the amount of quantum fluctuations 
in the system Jj, flj, [ij [iJI . This fact can be used as an 
operational definition of entanglement 0, Il3 ■ 

By the total amount of quantum fluctuations we mean 
the sum of variances of observables Oi which form an 
orthonormal basis of Lie algebra of the dynamic symme- 
try group G of quantum system [ll|. For example, in 
the case of two qubits (two two-level atoms, the dynamic 
symmetry group is G = SU(2) x SU{2), while the local 
observables are provided by the Pauli (spin projection) 
operators. 

The environment of Ref. [3 provided by the squeezed 
vacuum always manifests a high level of quantum fluc- 
tuations, and by means of interaction with the atomic 
system may transmit the quantum fluctuations to the 
atomic system, which is supposed to be initially unen- 
t angled. 

It seems to be interesting to consider the case when an 
entangled state is generated from initially unentangled 
state in the presence of non-fluctuating classical environ- 
ment. In this note, we show that a reasonable amount 
of entanglement can be achieved in a system of two-level 
atoms with dipole interaction in presence of a classical 
driving field. Although the classical driving field is de- 
void of quantum fluctuations, it helps to realize the po- 
tential ability of entanglement hidden in the structure 
of the cigcnstates of the system. Let us note that the 
stabilizing influence of the classical driving fleld on en- 
tanglement in a system of two-level atoms in a cavity 
with high quality factor has been discussed in Ref. |13l |. 



The paper is arranged as follows. In Sec. II, we discuss 
the model Hamiltonian and its properties. We examine 
the structure of eigenstates and show the presence of en- 
tanglement. We show that basis of the corresponding 
Hilbert space can be divided into the antisymmetric sin- 
glet part, represented by a certain maximum entangled 
state, and the triplet part. In Sec. Ill, we show that un- 
der certain physical assumptions the master equation can 
be solved in the triplet domain of the Hilbert space. Un- 
der the assumption of short interatomic distances (Lamb- 
Dicko limit), we find an exact steady state solution which 
give high enough amount of entanglement. We also ex- 
amine entanglement as a function of interatomic separa- 
tion and classical driving field. Sec. IV contains discus- 
sion of the obtained results and conclusions. 



II. PROPERTIES OF THE MODEL 
HAMILTONIAN 

The two-atom system under consideration can be spec- 
ified by the Hamiltonian 
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Here A^'^ = u)''^^ —ujc is the detuning between the atomic 
frequency lu^"^' and classical driving field frequency cue 
assumed to be small (w*^'^ '3> lA*^*)!), S*^*) represents the 
classical driving field at the location of the i-th atom, 
and the dipole-dipole coupling constant has the form [91 
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where /!''■* is the atomic dipole moment, ko 
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This Hamiltonian (1) is defined in the four-dimensional 
Hilbert space of two qubits 7^2,2 spanned by the vectors 
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where \af3) = \a) (E) \f3) and a = e,g denotes the excited 
and ground atomic states, respectively. This space 7^2,2 
represents in fact a product of the two two-dimensional 
"spin-i" spaces H1/2 spanned by the vectors |e) and g) 
each. This space can also be expanded as follows 



7^2,2 = 7^1/2 ® ^1/2 = Ho ® Hi, 
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where Hq is the singlet space, consisting of the anti- 
symmetric maximum entangled state 
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("spin-0" space). Using the terminology of the phase- 
transition theory [ij] , we can associate this state with the 
"aniferromagnetic" order in the system of two dipoles. 

In turn, the "spin-1" space Tii is spanned by the triplet 
of symmetric states 
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Formally, this basis (6) can be associated with the states 
of a "spin-1" object. Namely, the states | + 1), |0), and 
I — 1) correspond to the states with projections of "spin" 
equal to 1, 0, and —1, respectively. It should be noted 
that the state 10) of a qutrit always manifests maximum 
entanglement jlll[l3 |. 

Using the analogy with the phase transition theory, 
we can say that the unentangled states | -I- 1) and | = 1) 
correspond to the "ferromagnetic" (parallel) order of the 
dipoles, while the maximum entangled state |0) in (6) 
corresponds to the disordered "paramagnetic" state. 

Assume now that the atoms are identical when we may 
put A(i) = A(2) = A and r^^) = T^^) = r. In addition, 
we assume that the polarization of classical driving field 
is parallel to the interatomic axis, so that i?^^^ = i?^^) = 
E. It can be easily seen that under these conditions the 
antisymmetric singlet state (5) is an eigenstate of the 
Hamiltonian (1) with the eigenvalue ca ~ ~^'- 
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In the triplet sector, equation for eigenvalues takes the 
form 



e' - ne^ - (A^ + 4£;2)e + A^J] = 0. 
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In the absence of external field {E = 0), this equation 
has the roots e+i = A, eo = ri, and e_i = —A, and 
the set of the corresponding eigenstates is provided by 
the states (6). Thus, at fi > |A|, the "antiferromag- 
netic" maximum entangled state (5) is the ground state 
of the system in the absence of the external field. Oth- 
erwise, one of the states with "ferromagnetic" (parallel) 
order of atomic dipoles becomes the ground state. This 
is quite natural result, because the pseudo-spin Hamil- 
tonian (1) has the same structure as the model of the 
so-called plane rotator in transverse field, examined in 
the theory of phase transitions T4] . It is clear that sim- 
ilar behavior takes place at small values of the classical 
driving field. 

In the absence of detuning (A = 0) and an arbitrary 
value of the classical driving field, there is the following 
symmetric maximum entangled eigenstate 
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in the triplet sector, which corresponds to the eigenvalue 
Eq = 0, and the two eigenstates 
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with the eigenvalues 
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The concurrence of the states (9) is then given by the 
expression |l2| 
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where [tp±] denotes the matrix of coefhcients in the states 
(9). For the lowest state j-i/'-) with the eigenvalue e_ in 
(10), the concurrence (11) achieves the maximum value 
C = 1 at i? = and then monotonically decreases with 
the increase of E. Similar behavior takes place in the 
case of state \ip+) with the eigenvalue e+. 

Thus, the eigenstates of the system with the Hamilto- 
nian (1) can manifests entanglement and even maximum 
entanglement (at least at a certain choice of parameters) . 
Therefore, all one can assume is that the system under 
consideration has a good potential ability to evolve to- 
wards an entangled state under influence of a proper dis- 
sipative environment. 



III. IRREVERSIBLE EVOLUTION TOWARDS 

ENTANGLEMENT 



-afa'-I^p-pafcT^l^), (12) 



where the coefficients F^'^^ are responsible for the spon- 
taneous decay and have the form F^"^ = F^'^ and 
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We now show that under a certain choice of the initial 
state, the antiferromagnetic maximum entangled state 
(5) can be discarded from the evolution picture described 
by the master equation (12). Since (5) is the eigenstate 
of the Hamiltonian (1), it does not contribute into the 
unitary evolution, described by the first term in the right- 
hand side of (12). We also have 



^a«|A)=0. 



4=1 



Thus, it does not contribute into the Liouville term in the 
right-hand side of (12). This means that if the system 
is prepared initially in the maximum entangled antisym- 
metric state (5), it will stay in this state forever. If the 
initial state is chosen differently, the irreversible evolu- 
tion of the system should be considered in the triplet 
sector (6). 

This means that the density matrix in Eq. (12) takes 
the following block form 



The irreversible evolution assumes dissipation of en- 
ergy. In the system under consideration, it can be de- 
scribed by the losses due to spontaneous decay of the 
excited atomic states. We can specify the irreversible 
evolution in the system described by the Hamiltonian 
(1) with the assumptions made in the previous Section 
by the following master equation [Ig 
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Here the blocks correspond to the bases (6) and (5), 
respectively. The non-Hermitian effective Hamiltonian, 
corresponding to the evolution described by the master 
equation (12) with the density matrix (14), has the form 
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Here again F*^^) = F^^) = y. It is seen from the non- 
Hermitian Hamiltonian (15) that in the absence of the 



classical driving field, all states except \gg) are damped, 
so that the steady state entanglement at i? = is impos- 



sible, and the system evolves towards the unentangled 
state \gg). 

In other words, only the presence of classical driving 
field can stabilize steady state entanglement in the model 
system with the Hamiltonian (1). 

Because we are interested in the robust entanglement, 
let us consider the steady state solutions of the master 
equation (12) with the density matrix (14) in the triplet 
sector of the Hilbert space. 
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Consider first the Lamb-Dicke limit of short inter- 
atomic separation. Then, it follows from the definition 
of the decay rate (13) that 

p(12) ^p(21) ^Y. 

In this case, making the further assumption that A = 0, 
for the steady state density matrix in triplet sector we 
get the following explicit expression 
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Here N = Trp is the normalization factor and fi and E 
are established for the dimensionless parameters il/F and 
E/T, respectively. It follows from the numerical analy- 
sis that the detuning A has a weak influence upon the 
concurrence. 

To determine the settings, leading to the maximum 
possible amount of entanglement in the system under 
consideration, we choose Q, = tE"^, where r is a dimen- 
sionless constant to be determined upon the maximiza- 
tion of concurrence. This factor in the Lamb-Dicke limit 
can be represented as follows 
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where a — 1/137 is the fine structure constant, Q denotes 
atomic quality factor {Q = ujqT, where T is the lifetime 
of the excited atomic state), n is the mean number of 
photons per unit volume in classical driving field, and V 
denotes the volume of interaction between atom and field, 
so that nV gives the mean number of photons interacting 
with atom during the time T. 

The concurrence is defined as follows [la 



(18) 



C = max(Ai — A2 — A3, 0), 

where A denotes the spectrum of matrix R 
and p denotes the complex conjugation of (16) in the so- 
called "magic basis" [13 • The maximum entangled state 
provides C = 1, while the unentangled states give C = 0. 
One can see from Eq. (2) that at fixed r and in the 
Lamb-Dicke limit ko -f ^ 1 , both dimensionless parame- 
ters n/r, E/r ^ l. in this case, the density matrix (16) 
takes the form 
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To our surprise, the concurrence (18) in this limit turns 
out to be rational function of r 



extended by zero at r < 2. Thus, entanglement is impos- 
sible if T < 2. The maximum value of concurrence 
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0.43 
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is attained at 

Trr,ax^2 + 2VTS^9.21. 

The corresponding amount of entanglement \io^ is 
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0.285 ebit. 



Taking into account the form of the dimensionless pa- 
rameter T given by Eq. (17), we can examine the di- 
mensionless interatomic distance k^ ■ r, corresponding to 
the maximum entanglement provided by Tmax = 9.21, as 
a function of the number of photons fiV , which should 
obey the condition nV ^ 1 in the case of classical driving 
field. The dependence is given by Fig. 1. It is seen that in 
the case of mean number of photons nV ~ 10, the inter- 
atomic distance should be of the order of 10~^A (where 
A is the wavelength) to achieve the maximum possible 
amount of entanglement. Increase of the mean number 
of photons in the driving field, considered as a coherent 
state with jap ^ 1, leads to a decrease of interatomic 
distance, which is required to have maximum amount of 
entanglement. 

The dependence of concurrence from the interatomic 
distance comes from the coupling constant Vt (2). The 
results of numerical calculations for different interatomic 
separations and different values of the classical driving 
field are shown in Fig. 2. It is seen that deviation from 
the Lamb-Dicke limit decreases the concurrence. 

Thus, the classical driving field stabilizes entanglement 
in the system under consideration. 

We now turn to the discussion of the obtained results. 



IV. SUMMARY AND CONCLUSION 

We have examined the system of two identical two- level 
atoms, interacting with each other by means of dipole 
forces. The dissipation of energy in the system is pro- 
vided by the spontaneous decay of the excited atomic 
states. The compensation of losses is provided by the 
presence of classical driving field. 

The consideration of the eigenstates of the Hamil- 
tonian (1), describing the system under consideration, 
shows that the system carries a potential ability of cre- 
ation of entanglement, hidden in the structure of the 
eigenstates. 

Instead of conventional picture of two-qubit system, 
we prefer to specify the system by the triplet of states 
of an effective qutrit and by an "antiferromagnetic" , an- 
tisymmetric maximum entangled state, corresponding to 
a "spinless" quasiparticle. 

It is shown that in the absence of the classical driving 
field, system evolves towards an unentangled state (both 
atomic dipoles are in the ground state). The presence of 
the classical driving field stabilizes the entanglement. 

In the Lamb-Dicke limit of a point-like system, we 
obtained an exact solution for the steady state den- 
sity matrix, that manifests high amount of entanglement 
(the concurrence Cmax = 0.43 and the entanglement 
£max = 0.285 ebit). This amount is much higher than in 
a number of recent proposals. In particular, it is higher 
than that in the case when the squeezed vacuum is used 
for stabilization of entanglement instead of the classical 
driving field Q- It should be stressed that this amount 
of concurrence corresponds to a reasonable magnitude of 
the parameter _E, describing interaction of classical driv- 
ing field with atoms. 

Since the classical driving field is devoid of quantum 
fluctuations, the only source of the high-level quantum 
fluctuations peculiar for entanglement can be associated 
with the structure of the eigenstates of the Hamiltonian 
(1), describing the system. Thus, the classical driving 
field plays a double role. Namely, it compensates the 
losses of energy and stabilizes the entanglement carried 
by the system. 

Such a picture suggests an analogy with the phase tran- 
sitions in antiferromagnets, whose ground state can be 
considered as a collective entangled state. This ground 
state of such a system is realized at low temperatures. 
In other words, it is stabilized by the presence of a low- 
temperature environment (thermostat). 

In our case, the system under consideration has very 
few degrees of freedom and therefore cannot manifest a 



thermodynamic phase transition. Nevertheless, the clas- 
sical driving field together with the spontaneous decay 
plays a role similar to that of the thermostat. 

In the investigation of steady state entanglement in 
the system, we always considered the triplet sector of 
the Hilbert space. From the physical point of view, this 
means that the initial states are symmetric. For example, 
the most natural choice is provided by the unentangled 
initial state \gg) with both atoms in the ground state. 
In this case, a certain balance between the spontaneous 
decay and influence of the classical driving field lead to 
creation of the steady state entanglement in the system. 

Let us now discuss the effect caused by the inclusion of 
the maximum entangled antisymmetric state (5) into the 
picture. This means that we choose the density matrix 
of the form 

p\A){A\+qps, 

in the whole four-dimensional Hilbert space. Here ps de- 
notes the block, coming from the triplet sector and p, q 
are the weights of the corresponding contributions (p -I- 
q = 1). Since the state \A) belongs to the "magic basis" 
|l5| , the concurrence of the steady state is given by for- 
mula similar to (17) applied to spectrum q\i, q\2, qX3,P- 
As far as p < qXi, that is p < Ai/(1 -|- Ai), the pres- 
ence of the antisymmetric component (5) decreases the 
concurrence. 

Thus, to achieve a reasonable amount of robust entan- 
glement, we have to prepare the system initially either in 
the antisymmetric maximum entangled state (5) or in a 
state from the symmetric sector (6) of the Hilbert space. 

In the above consideration, we always assumed that 
atoms are identical. It seems to be interesting to extend 
our consideration on the case of non- identical atoms. In 
view of the result of Ref . Q , we can expect that this may 
lead to a significant increase of amount of entanglement. 

We also restricted our consideration to the case of po- 
larization of the classical driving field parallel to the in- 
teratomic axis. The alternative choice of the polarization 
perpendicular to the interatomic axis can lead to a strong 
changes of picture as well. First of all, the change of po- 
larization changes the form of the coupling constant (2). 
Then, it causes the consideration of the different values 
of the classical driving field in the atomic locations. Fi- 
nally, in this case, the antisymmetric state (5) cannot be 
discarded at all, so that the irreversible evolution should 
be examined in the whole four-dimensional Hilbert space. 

The detailed analysis of the above mentioned two ex- 
tensions of the model deserves special consideration. 
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FIG. 1: Lamb-Dicke parameter (fco ■ f) versus the number of 
photons in the case of dimensionless parameter (17), provid- 
ing the maximal amount of entanglement for typical atomic 
quality factor Q = 10~^. 



FIG. 2: Numerical dependence of concurrence on the inter- 
atomic distance and classical driving field. The dimensionless 
quantities fl/T and E/T are used here. 



